Scaling and the center of band anomaly in a one-dimensional Anderson model with 

diagonal disorder 
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We resolve the problem of the violation of single parameter scaling at the zero energy of the 
Anderson tight-binding model with diagonal disorder. It follows from the symmetry properties of 
the tight-binding Hamiltonian that this spectral point is in fact a boundary between two adjacent 
bands. The states in the vicinity of this energy behave similarly to states at other band boundaries, 
which are known to violate single parameter scaling. 



INTRODUCTION 

Despite a long history of research in the field of lo- 
calization, even such a standard problem as the one- 
dimensional Anderson tight-binding model with nearest 
neighbor hopping and diagonal disorder can still hold sur- 
prises. This model is described by the equation of motion 



ip n+1 + ^> n _i + (e„ - E)ip v 



0, 



(1) 



where, tp n represents a wave function at the n-th site, 
e n is a random site energy distributed uniformly between 
— W and W, and E is the energy. Recently, Schomerus 
and Titov Q demonstrated that single parameter scaling 
(SPS), which is a cornerstone of the conceptual frame- 
work for dealing with transport properties of disordered 
systems 0, is violated in this classical model of local- 
ization in the vicinity of the energy E = 0. While the 
violation of SPS in the vicinities of the energies E = ±2 
(which are band boundaries of the system with W = 0, 
i.e., all e„ = 0) is well established 0, the absence of 
SPS at the band center E = looks surprising. It is yet 
another manifestation of the anomalous nature of this 
spectral point. Whereas earlier studies of the Anderson 
model revealed an anomalous behavior of the localization 

length at £ = HHH0], the violation of SPS in this 
spectral region remained unnoticed. The mathematical 

roots of the anomalous behavior of the localization length 
were established in Ref. 0. At the same time the physical 
understanding of the nature of the anomalous properties 
of E = is still absent. Here we suggest an approach to 
this spectral region based on symmetries of the Hamil- 
tonian Eq. (JIJ j which clarifies the physical nature of this 
anomaly, and explains scaling properties of the Anderson 
model in this spectral region. 

It is convenient to describe the statistics of conduc- 
tance, g, of a disordered chain of length L in terms of the 
Lyapunov exponent (LE), 7, defined as |j| 



7 =iln( ! 



(2) 



distribution can be expressed through the mean value of 
LE, 7 = (7). Such a relationship between 7 and the 
variance, a 2 = (7 2 ) — (7) 2 , was first conjectured in Ref.[|| 



T = (T 2 i/7 = 1. 



(3) 



SPS means that the probability distribution of 7 depends 
on a single parameter, so that, e.g., any moment of this 



The authors of Ref. |jj demonstrated that at E — 0, pa- 
rameter t of Eq. 10 deviates from the SPS value and is 
equal to 1.047. This deviation from Eq. ©, which was 
overlooked by previous studies probably due to its small 
size, is of principle importance. It seemingly contradicts 
a well established criterion for SPS introduced in Ref. 3. 

In this Letter we show that this is only an appar- 
ent contradiction, and explain the violation of SPS at 
E = within the general picture developed in Ref. 0. 
The anomaly at E — turns out to result from a hidden 
symmetry of the tight-binding model, which splits the 
conduction band into two adjacent bands with a bound- 
ary at E = 0. This splitting does not affect thermody- 
namic properties of the model but is crucial to its scaling 
behavior. 

Before presenting our resolution of the mysterious be- 
havior in the vicinity of E = 0, let us recall the histor 
of the criterion for SPS. Originally, Anderson, et al. 
suggested that SPS holds when the localization length 
of the system, l\ oc — 7" 1 , exceeds what they called the 
phase randomization length: ij oc > l p h. The stationary 
distribution of the phases of the reflection and transmis- 
sion amplitudes in disordered one-dimensional systems 
was assumed to become uniform as soon as the system 
length L exceeds l v n ■ Many authors derived Eq. © from 
this phase randomization hypothesis for several different 
models. The phase randomization itself was rigorously 
proven only for few particular models: banded matri- 
ces |9|, the Anderson tight-binding model with diagonal 
disorder 0, and a continuous model with a white-noise 
random potential for certain spectral regions. At the 
same time, it was shown that for some values of energies, 
such as E = or E = ±2 for the model described by 
Eq. (JTJ the stationary distribution of the phases is not 
uniform 0, S 0, [3 El 

Although in the vicinities of these special energies SPS 
is indeed violated 0, 0], the criterion for SPS based 
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on the phase randomization length does not seem to de- 
scribe the situation adequately. Indeed, the absence of 
the phase randomization in all situations discussed above 
does not mean that the "relaxation" length of the phase 
distribution diverges and exceeds the localization length 
when one approaches the respective energies. What ac- 
tually happens is that the uniform contribution to the 
stationary phase distribution vanishes in the vicinity of 
these special energies. As a result, the nonuniform part, 
which always exists but is usually small, becomes domi- 
nant [M 0, 0] - Thus, the phase randomization length, as 
it was introduced in Ref. Q, is ill-defined. 

Thus, appealing to the phase distribution in connection 
with SPS simply substitutes the problem of formulating a 
criterion for SPS by the problem of finding a criterion for 
the uniform distribution of phases. Previously proposed 
criteria for the latter are specific for particular models. 
For instance, in Ref. Q] for the Anderson model within 
the white noise approximation, the criteria for the phase 
randomization for the band edges and the band center 
are set as \E - 2| > 3W 2/3 and \E\ > 10D, respectively 
(D is a measure of the disorder strength). These defini- 
tions do not look satisfactory to us. We believe that the 
criterion should be universal and formulated in terms of 
some fundamental macroscopic quantities. 

Such a universal criterion for SPS similar in form to the 
Anderson's, k oc > l p h, was recently suggested in Ref. 0: 

K = ll oc /l s > 1. (4) 

The new length scale, l s , has nothing to do with l p h and 
can be expressed through the integral density of states 
normalized by the total number of states in the band, 
N(E), and a distance between the neighboring sites, a: 

l s = a/sm[irN(E)]. (5) 

This criterion was extracted from the exact calculation 
of the variance of LE for the Anderson model with the 
Cauchy distribution of the site energies, e„ (the Lloyd 
model) . Eq. © was derived in Ref. y without any rela- 
tion to the phase distribution [l^. For in-band states, 
N(E) ~ 1/2 and l s is microscopic and therefore is not sig- 
nificant. For the states close to the band edges and fluc- 
tuation states in the former band gaps, where N(E) <C 1 
or 1 — N(E) <C 1, l s becomes macroscopic and can exceed 
k oc . In the latter case, l s has a clear physical meaning 
as an average distance between localization centers re- 
sponsible for the states with energies between E and the 
closest fluctuation boundary of the spectrum. 

The expression of l s in terms of the density of states 
N(E) allowed authors of Refs. to conjecture a gener- 
alization of this definition to other models. It was sug- 
gested that N(E) in Eq. © must be understood as the 
total number of states (per unit length) between E and 
a closest genuine spectral boundary normalized by the 
total number of states. This definition of l s can be used 



also in the case of systems with multiple bands, provided 
that the disorder does not fill the band gaps of the origi- 
nal systems completely, and there exist genuine spectral 
boundaries inside each of the band gaps. Under this con- 
dition, N(E) refers to the number of states associated 
with one particular band of the spectrum, and varies be- 
tween and 1 when E spans the states of the band. 

In light of this criterion, one can understand the vi- 
olation of SPS in the band gaps of the initial spec- 
trum observed numerically in Ref. Applicability of 
the / s -based criterion was examined for several other 
models such as the Anderson model with the box Q 
and dichotomic 0] distributions of the site energies, 
the model of a scalar wave propagating through a one- 
dimensional disordered periodic-on-average superlattice, 
and even for the model of a wave propagating in an ab- 
sorbing medium 0] . In all these cases, the criterion per- 
fectly described the transition between the SPS behavior 
of in-band states, and non-SPS behavior of the band edge 
and fluctuation states. Moreover, as was shown recently 
in Ref. lid, the length l s plays an even more important 
role than simply determining a boundary between SPS 
and non-SPS spectral regions: the second and the third 
moments of the distribution of the Lyapunov exponent 
in the non-SPS region can be parameterized by a single 
parameter k provided that L exceeds l s |l8j . 

THE BAND-CENTER ANOMALY IN THE 
ANDERSON MODEL 

The success of the scaling approach based on the pa- 
rameter k, Eqs. (@J) and in predicting violations of 
SPS at the band edges of different models, and in describ- 
ing the probability distribution of LE in non-SPS spectral 
regions 0, 0, motivated us to look more carefully at 
the reasons for the apparent failure of this approach at 
the center of the band of the Anderson model. 

We start our analysis with a tight-binding model with- 
out disorder, e n = 0. The solution of Eq. at e„ = 
is plane waves, ip n cx exp(±i/cno), with fc satisfying the 
dispersion equation 

£ = 2cosfca; < ka < tt, (6) 

According to the traditional widely accepted point of 
view, Eq. © describes a single energy band —2 < E < 2, 
However, this model contains more than the simple 
picture reveals. Indeed, the equations of motion, Eq. 
at e n = acquire an additional symmetry: Operation 

D n = V>„ -> (-l)"Vn (?) 

transforms a state ip n (E) into the state ip n (—E): 
tp n (E) = (— l) n ijj n {— E). This operation forms a group, 
and the states of this model can be classified according 
to its irreducible representations. 
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In order to realize these representations, we introduce 
two field invariants with respect to D n : 



ui = ^2 



Vl = 1p2l+l- 



(8) 



Solving equations of motion for these fields we find that 
the spectrum is separated into two branches described by 
the following dispersion equations 



E = ±cos(ka/2). 



(9) 



Two bands of the spectrum described by Eq. © have 
a common boundary at E = 0, and, therefore, the spec- 
trum of the system appears as consisting of a single band, 
as it is usually assumed. However, the classification of the 
states according to the irreducible representations of D n 
reveals that E = should be considered as a band bound- 
ary rather than as a band center. The new description of 
the band structure can be viewed as a transition to the 
reduced zone representation, where the wave number of 
the original dispersion equation Eq. © is restricted to 
the interval < ka < n/2. This operation is equivalent 
to doubling of the elementary cell of the original periodic 
chain. This can be justified by noticing that two adjacent 
sites of the chain are not equivalent with respect to the 
operation D n [Eq. J7J]. 

The considerations presented above would have been 
merely an exercise, if the vicinity of E — did not behave 
anomalously in the presence of disorder. Of course, any 
disorder mixes states from E < and E > bands. 
However, since violation of SPS occurs in the vicinity of 
a spectral boundary of an unperturbed system jjj, one 
can expect that E = behaves similarly to E = ±2 
boundaries. Actually a similarity between properties of 
E = and other band boundaries was noted already in 
Ref. 6, where it was shown that the calculation of the LE 
in the vicinities of all three points E = ±2, and E = 
required almost identical mathematical approach. 

The deviation of the parameter t(k), Eq. from 
unity observed in Ref. |l| can be interpreted as a mani- 
festation of the band-edge nature of the spectral point 
E = 0. This deviation is, however, different from the be- 
havior of t(k) in the vicinity of E = ±2. Near E = ±2, 
t(k) demonstrates a small overshoot above unity at the 
band side of these boundaries, and decreases quickly at 
the band gap sides Q. Since two conduction bands are 
adjacent at E — 0, there is no band gap related decrease 
of t(k). However, the band-edge nature of E = still 
manifests itself in the form of the overshoot, discovered 
in Ref. 1. This overshoot appears at both sides of E = 0. 

The fundamental question is whether the deviation 
from SPS near E = occurs in accordance with the 
criterion Eq. The answer is "y es '\ however, the def- 
inition of the length Z s , Eq. ©, should be modified to 
reflect the fact that the states are now distributed be- 
tween two bands. Accordingly, when calculating l s for 
each of the bands, we have to keep in mind that N(E) 



1.2 



1.1 



2 1.0 



0.9 



0.8 



0.7 



0.6 



x(k) near E = 
t(k) near band edge 



S 0.9- 



0.7 1 



-1.0 -0.8 -0.6 -0.4 -0.2 



FIG. 1: Dependence of the scaling parameter r on re in the 
vicinity of the fluctuating band edge \E\ ~ 2 and near E = 
for different strengths of disorder: W changes from 0.15 to 
0.4. Curves corresponding to different strengths are indis- 
tinguishable. The inset shows r(re) for the Anderson model 
with next to the nearest neighbors interaction, Eq. 110L in 
the vicinity of the energy E — — 2a for a = 0.3. 



in Eq. (J5J is a number of states between E = and a 
given energy E normalized by the total number of states 
in this band. 

The results of the numerical calculations based on the 
new definition of l s reveal first of all (Fig. ^| that the cri- 
terion k > 1 and SPS are indeed violated simultaneously 
near E — as well as near E = — 2. What we find even 
more important, however, is that the results of the cal- 
culations for different energies, strengths of disorder, and 
lengths of the system, all collapse to a single curve, when 
presented in the form r(re). Thus n = li oc /l S i after the 
proper definition of l s , remains a natural scaling variable 
for the entire spectrum of the system. 

This means that separating negative and positive en- 
ergies in two bands is physically significant. Of course, 
functions r(re) for the vicinities of E — —2 and E = 0, 
do not coincide. This does not cause problems for our 
scaling description as long as these two regions are well 
separated by the spectral interval where r = 1. 

Additional arguments in favor of our interpretation of 
the E = anomaly can be obtained by modifying the 
initial Anderson model. For example, one can include 
next to the nearest neighbors interaction, and consider a 
model determined by the equation 

ip n+ l + ip n -i + a(V>„+2 + ipn-2) + (e« - E)ip n = 0. (10) 

The spectrum of this model, in the absence of disor- 
der, consists of two branches, only one of which cor- 
responds to propagating states for a < 1/4. This 
branch is invariant with respect to a similar symmetry 
(— l) n ijj rL {E) = ip n (E'). The stationary point of this 
transformation E' = E = —2a is a boundary between 
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FIG. 2: Effect of the periodic potential on the scaling param- 
eter r near E = for W — 0.4 Squares, circles and triangles 
correspond to Uo = 0.002, 0.02 and 0.2, respectively. 



nature of this energy, and the center-band anomaly in 
the behavior of the second moment is not observed. 

Concluding, we resolved the mystery of the anomalous 
properties in the vicinity of E = in the Anderson model 
with diagonal disorder. We demonstrated that due to an 
additional symmetry of the respective Hamiltonian, this 
spectral point must be considered as a band boundary 
between two adjacent bands rather than the center of a 
single band. The relevance of the criterion for SPS sug- 
gested in Refs. 3, is thus reestablished. Our approach 
allows one to describe statistical properties of conduc- 
tance in the vicinity of E = within the framework of 
the scaling approach of Ref. 

The authors thank Steve Schwarz for reading and com- 
menting on the manuscript. The work at Queens College 
was supported by AFOSR under Contract No. F49620- 
02-1-0305, and partially by PSC-CUNY grants. BA 
gratefully acknowledges partial support by the EPSRC 
grant GR/R95432. 



two adjacent bands. Numerical calculations confirmed 
that indeed there is an anomaly in the dependence of 
the scaling parameter t(k) at —2a similar to the one 
discussed above (see inset in Fig. P). 

To justify the physical significance of our treatment of 
E = as a band edge further we added a periodic poten- 
tial U n = Uo(—l) n to random site energies e„. This po- 
tential breaks the symmetry of the model, and creates in 
the vicinity of E = a real band gap with a width propor- 
tional to Uq. The idea is to look at the sensitivity of the 
anomalous behavior in the vicinity of E = to the pres- 
ence of the real band gap. This gap should decrease to- 
gether with the magnitude of the symmetry-breaking po- 
tential. Will the features associated with the gap survive 
the vanishing of the gap width? The results of the nu- 
merical simulations answering this question are presented 
in Fig. |2 One can see that the anomalies in t(k) remain 
stable with respect to the symmetry-breaking perturba- 
tion. In the limit Uq — > 0, t(k) has the same behavior as 
for small but finite Uq. This observation in conjunction 
with the scaling arguments presented above gives strong 
support to our treatment of E = as a band edge. 

Finally, there is one more question that should be ad- 
dressed. Choosing an appropriate periodic potential one 



can create a band gap centered not only at E g 



but 



anywhere in the spectrum. Why, then, is SPS violated 
only near E = 0? It follows from our numerical studies 
that when the periodic potential is weak and E gap ^= 0, 
the 7"(/c)-function always approaches the band edge with- 
out any overshoot above its SPS value, r = 1; it simply 
decreases inside the gap. As a result, the normal SPS-like 
behavior is restored once the periodic potential vanishes. 
In the Lloyd model, the similar behavior is characteris- 
tic for all band boundaries, E = ±2, 0. Therefore, the 
absence of a real band gap at E = masks the actual 
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